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The empirical-phenomenological quasi-deuteron photo-fission description is theoretically justified 
within the semiclassical, intermediate statistics model. The transmutational fermion (nucleon) - 
boson (quasi-deuteron) potential plays an essential role in the present context and is expressed in 
terms of thermodynamical and of microscopical quantities, analogous to those commonly used in 
the superfluid nuclear model. 
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I. INTRODUCTION 



Over the past several years, many efforts have been directed to develop methods for the physical description of 
complex many-body systems. The implementation of quantum approaches like the Hartree-Fock, the Random Phase 
Approximation and the Time Dependent Hatree-Fock Method have allowed the explanation of a great variety of 
experimental results and particularly of the properties of the nuclear ground state [1,2]. 

In the literature, the nuclear dynamics is described by different theories. In the most cases, approaches based on 
semiclassical one-body equations (Fokker-Plank, Boltzmann, Boltzmann-Langevin equations) are use. These methods 
f* — | can be regarded as semiclassical analogs of the mentioned quantum approaches [3,4] and include the Pauli exclusion 
and the inclusion principle [5] to represent correctly systems of fermions, bosons or anyons. 

In several problems in condensed matter and in nuclear physics, the correlation effects between pairs of fermions 
are quite relevant in the interpretation of experimental results. Similarly, the interactions among bosons are relevant 
to the superfluid nuclear model, the interacting boson model and the mean field boson approximation and allow the 
Q\ • explanation of many collective nuclear properties. In these cases, the interaction among the valence nucleons outside 
the core produces pairs of correlated nucleons that can be approximated with bosonic particles (like, for instance, the 
> quasi-deuterons). The pairing of fermions leads a certain amount of nucleons to a transition from one statistics to 
another (e.g. Bose-Einstein). This transition can be modeled with a transmutational potential (TP) and implies that 
a nucleon can be considered, in particular excited states, as an ensemble of fermions and bosons or as an ensemble of 
particles of intermediate statistics. 

Semiclassical kinetic approaches to intermediate statistics [5-7] are based on the classical Fokker-Planck equation 
and describe the time evolution of the distribution function n(t, v). The Fokker-Planck equation is corrected with 
the insertion of the factor 1 + nn(t,v) into the transition rates, where k varies between —1 (Fermi-Dirac distribution 



(FD)) and +1 (Bose-Einstein distribution (BE)), (n = is the case of the Maxwell-Boltzmann distribution). The 
effect on the distribution function of the exclusion (k < 0) and of the inclusion (k > 0) principle, can be evaluated 
with this approach. 

This approach is not unique, different non linear expressions can also be considered [7]. However, more involved 
expressions are not required for the discussion presented in this work. 

In the present paper, we call statistical transmutation the variation of the particle state from its characteristic value 
k to a new value k' (this definition includes, in particular, the transmutation from k = —1 to n' = +1). 

To describe the importance of the TP, we fix our attention to the photofission of heavy nuclei in the energy range 
of excitation between 40 MeV and 100 MeV (usually called quasi-deuteron (QD) energy region [8,9]). 
In photonuclear reactions, it is well known that, above the GDR and below the pion photoproduction threshold, 
the photoabsorption involves two nucleons to satisfy the kinematical conservation laws. This circumstance helps the 
formation of strongly correlated pairs of nucleons (quasi-deuteron states). The photoabsorption cross section is pro- 
portional to the deuteron photodisintcgration cross section. 
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Recently, we have developed an extension of the QD model, to describe the total photoabsortion and the photofission 
of heavy nuclei [10-11] (an extension to higher excitation energies, pion photoproduction region, of the QD photofis- 
sion model has been recently developed by Arruda-Neto and collaborators [12,13]). The agreement of the model 
with the experimental results is quite good for a wide variety of nuclei (Bi, Th, U) (see also Ref.s [14-16]). Our 
phenomcnological-empirical model provides us useful quantities, such as the number of QD participating to photofis- 
sion at different excitation energies. Its use, together with the semiclassical kinetic approach, provides a theoretical 
basis for the phcnomenological model itself. The evaluation of the TP in analogy with the superfluid model, allows 
the introduction of an analytical expression for the excited QD nuclear state as a mixture of bosons and fermions. 
The TP can be given an interpretation in terms of the occupation operators. An interpretation of the TP in terms of 
thermodynamical quantities is also given. 

In Sect. II, we describe the equations of the transmutational kinetic model and provide the free nucleon fraction 
and the probability of participation of one nucleon to a QD pair. 

In Sect. Ill, we calculate the transmutational potential using the results of a phenomenological approach to photofis- 
sion, valid in the QD energy region and expressed in terms of the nuclear latent heat. 

In Sect. IV, we give an interpretation of the transmutational potential using an analogy with the superfluid model of 
the nucleus. Conclusion are outlined in Sect. V. 



II. NUCLEON-QUASI DEUTERON TRANSMUTATION MODEL 

The protons and the neutrons are distinguished by their third isospin component t 3 — r (f = — r). The nucleon 
mass is indicated by m N = 939 MeV and the occupational number in the velocity space by nT(i,v). The QD 
particles, neglecting their binding energy, have mass m D = 2m N , L = 0, S = 1, T = 0. Their occupational number 
is n D (t,v). 

By using an intermediate statistics with a linear enhancement or inhibition factor to calculate the QD photofission 
cross section, we can assign to the nucleons a value of the parameter k, intermediate between —1 and 1 and different 
from zero. When the excitation energy increases, the parameter n varies because the QD particles are more numerous 
as the nucleus is more excited. Therefore, the bosonic component of each nucleon increases and one can study the 
TP from the statistical transmutation from k to k . We are interested in the transmutation nucleon-QD or viceversa, 
therefore we must fix, in this context, the states n and k as the fermion and boson state, respectively. 
An alternative, but equivalent, procedure is to assume that the system (nucleus) is composed by a certain number 
of fermions (nucleons) and by a certain number of bosons (QD), depending on the excitation energy. Therefore, we 
write two different kinetic equations for the nucleons, with k = — 1, and for the QD, with k = +1. In this way it can 
be clearly understood the role of the fermion-boson TP. 

We consider the energy as a continuous variable, the two functions n T N (t, v) and n D (t, v) obey the following system 
of coupled equations [6]: 

^^ + vj;(t,v )+ . ? ;;„ D (t,v) = o , (i) 



+ Vj D (t, v) + jZ 2N (t, v) = . (2) 

These equations can be viewed as the continuity equations of a nucleon and of a quasi-deuteron gas in the velocity 
space. When the energy is a discrete variable, we can substitute Eq.s (1) and (2) with a system of coupled master 
equations. 

In the velocity space, the currents (t, v) and j D (t, v) have the following expressions: 

j; (t, v) = - [j; (t, v) + VDl (t, v)] nl (t, v) [1 - nl (t, v)] - IT (t, v)V< (t, v) , (3) 



j D (t, v) = - [J D (t, v) + VD D (t, v)] n D (t, v) [1 + n D (t, v)] - D N (t, v)Vn D (t, v) , (4) 

where 3 NjD (t, v) and D N D (t,v) are, respectively, the drift and diffusion coefficients for the nucleon and the QD. 
The net transmutational current J (J^I^n ) takes into account the formation of a quasi-deuteron (two nucleons) 
from two nucleons (a deuteron) with different third isospin component; its expression, according to the exclusion- 
inclusion principle (EIP), is the following [6]: 
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f 2 Lo (*. v ) = r 2N^o (*) 9 N (*, v) n T N (t, v) n f N (t, v) [1 + n D (t, v)] 

-r D ^ N (t)g D (t,v)n D (t,v)[l-n T N (t,v)][l-nl(t,v)} , (5) 

where r 2N ^ D (t) is the transmutation rate of two nucleons in a quasi-deuteron, r D ^ 2N (t) is the transmutation rate of 
a quasi-deuteron into two nucleons, g N (t, v) and g D (i, v) arc functions taking into account the interaction among the 
identical particles belonging to the bound system we are considering. 
From the definition of the net transmutational current we have: 

3lLo(^) = -.f D l 2 J^) (6) 

Equations (1-6) define univocally the diffusion process of neutrons, protons and quasi-deuterons in the velocity space 
and the formation and disgregation of quasi-deuterons inside the nucleus. 

The statistical distribution n N D (v) can be obtained, in stationary conditions (i — > oo), when the currents (t, v) 
and j D (t,v) vanish. 

If D n,d(v) = D N D {v) and J N _ D (v) =\J N D {v)/v wc obtain: 

n N,nM = Vm~P " VI ' ( ? ) 



where E N D = \m N D v 2 + V N D (v) and 

dV N:D (v) 



8v PD ND (v) 



(8) 



k n = —1, k d = 1, (3 = 1/kT and T is the temperature of the system, fj, N D is the chemical potential of the nucleons 
and the QD particles. In Eq. (8), dV N D (v)/dv represents the difference between the force acting among the particles 
and the Brownian force m N D v. 

We note that, in the case of Brownian particles, the potential V ND (v) is a constant and n N D (v) reproduces the 

standard FD and BE statistical distribution. 

Wc define the transmutational potential r] 2N ^ D as 

!^=ex P [/^ D ] (9) 

and introduce the potential h(v) 



9 D (v) 



' f,) = e X p[(3h(v)} . (10) 



The condition that the net transmutational current, given by Eq.(5), must satisfy the equation j 2NMD (oo,v) = 0, 
imposes the two following relations: 

V2N^D = M D - 2 M» > (11) 



h(v) = V^(v) + V^v)-V D (v) . (12) 

Equation (11), at any excitation energy E 11 provides a relationship between the chemical potentials of the two different 
types of particles and the TP. Eq. (12) allows us to obtain the potential h(v) in terms of the interaction potentials 
(constant in the case of brownian particles) . 
Let us define the two quantities: 

t».» = > ( 13 ) 

which are, respectively, the free nucleon fraction and the quasi-deuteron fraction, A is the mass number of the nucleus. 
In other words: £ N and 2£ D represent the probability that a nucleon is a free nucleon or a part of a QD. £ N and £ D 
satisfy the normalization condition: 

C+2^ = l • (14) 
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The fraction £ N D of nucleons and of QD is given by: 



= f N ,o ~ p (^) 3 / ^nMd'v , (15) 

where p is the nuclear density and f N D is the spin-isospin degeneration factor (f N = 4 for the nucleons; f D = 3 for 
the QD, having total spin 5=1 and total isospin T = 0), V is the nuclear volume, h the Planck constant. 
Note that we have considered pairs with L = 0, S = 1 and T = 0, that are the quantum numbers of a free deuteron. 
Recently, it has been considered the presence of QD pairs with T = 1 to obtain agreement with the experimental 
results on photoabsorption [17-19]. In our case, if we considered a QD component with T = 1 and S = 1, we should 
have taken the degeneration factor f D = 6. In that case, the value of the QD chemical potential would be reduced, 
because the total number of particles is fixed by the Eq.(15) and also the value of the transmutational potential r) 2N ^ D 
would be reduced because of Eq.(ll). The insertion of this component is not relevant to the results obtained in this 
work. 

If, for a given nucleus, £ D and £ N are known, as a function of the excitation energy, inverting Eq.(15), one can 
deduce the chemical potentials fi N e \i D and, by means of equation (11), derive the transmutational potential rj 2N ^ D 
as a function of £L . 



III. CALCULATION OF THE TRANSMUTATIONAL POTENTIAL 



If the energy of the incident photon is fully converted into excitation energy, the nuclear temperature T can be 
calculated from the well known relation: 



£L = — T 2 MeV 



(16) 



In our case this relation is not rigorously exact, as shown by Montecarlo calculations in an intranuclear cascade model 
[20,21]. Indeed, not all the photon energy contributes to increase the nuclear temperature. In the QD energy region, 
the excitation energy is about 15% smaller than the photon energy [21]. We have verified that our results are not 
modified if we consider this reduction (for a discussion on the validity of Eq.(16) see Rcf.[22]). 

The QD and nucleon fractions have been determined using our QD model of photofission [10]. Let us recall that 
the following quantities. First, 



NZ 

C(N,Z) = 7.72 — 



(17) 



is the effective number of QD pairs used to write the photabsorption cross section as the product of the QD effective 
number times the photodisintcgration cross section of the deuteron a D . Second, 



F 1 {E 1 



exp 



D 



(18) 



is the probability that one of the C (N, Z) pairs takes part in the photoabsorption reaction (this factor is due to the 
Pauli blocking and the quantity D/2 is the average energy to excite each nucleon of a pair above the Fermi level 
(D = 60 MeV) ). Third, 



F 2 {E 1 ) =cxp - 



D - T(E 7 ) 

E<y 



1 — exp I — 



D + TjE^ 



(19) 



is the probability that one of the C(N, Z) pairs participates in a photoabsorption evolving toward photofission (the 
product C(N, Z)F 2 a D is the photofission cross section); T(E 7 ) is a polynomial in E 1 (we have introduced it to reduce, 
in the first factor of F 2 , the photofission probability at low excitation energies, where the photoabsorption without 
fission is the most probable process, and to take into account, by means of the second factor, the photofission reduction 
at high energies). In the QD energy region, the contribution of the photabsorption to the continuum is greater than 
other contributions. T(E 7 ) represents the window of energy states, around the Fermi level, allowed to be occupied by 
the excited nucleons and pairs to arrange a system of particles evolving toward fission. 
From Eq.s (13) and (14) the number of QD particles, correctly normalized is: 
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1 

H N (E 1 ) = l-F 2 (E J ) . (21) 



tn(E 7 ) = ME,) , (20) 



We have considered three nuclei with very different photo-fission features: 209 Si, 232 Th, 238 U. 
In the Table the quantity F^E^) for the three different nuclei is reported at different excitation energies. The 
numerical values are obtained from the functions T(_E 7 ) given in Rcf.[10]. 

We can calculate the QD and the nucleon fractions, at any excitation energy, with Eq.s (20,21). Using Eq.(15), we 
deduce the QD and nucleon chemical potentials and finally, using relation (11), the TP. 
The nuclear density has been assumed equal to the infinite nuclear matter, i.e. 0.17 nucleon fm~ 3 . 
The nuclear interaction of Eq.(12) is described by a square well with V N = V D = V = —48 MeV. At the low 
excitation temperature T = 1.1 1.9 MeV, the chemical potentials fj, D w V , for 238 U and 232 Th and varies between 
-55 MeV and -51 MeV for 209 Bi. 

The quantity ?7 2JV ^ D (-E 7 ) is an increasing function of E 7 , as shown in the Figure. As the quantity of n 2N ^ D 
approaches a positive value, the number of quasi deuterons participating to the photofission becomes greater than the 
number of single nucleons. Infact, for 238 U, n 2N ^ D (90 MeV) — and at this energy F 2 (90 MeV) w 0.5, which equals 
the value reported in the Table. Only in the case of 209 Bi, n 2N ^ D (Ej) is nearly constant and negative. Let us note 
that, when n 2N ^ D — > — oo, the full system is composed by fermions (see Eq. (9)), on the other hand, if rj 2N ^ D — > +oo 
the full system is composed by bosons. 

In ref. [10,11] we have shown that the photofission probability Pf — Pf{E~ f ) is the solution of the following equation 



r * ( 22 ) 
i f E 7 

where f(^/~E 1 ) is a polynomial in \[E^. 

When f{\fE 1 ) is nearly constant (as in the case of 209 Bi), the solution of Eq.(22) is given by the simple expression 
Pf = a + bj\[E 1 which is the result of the statistical model. Usually f(\f~E 1 ) is not constant and the expression of 
Pf is more involved. 

Eq.(22) can be seen as a generalized Clausius-Clapeyron equation where f{\/E 1 ) plays, in the photofission, the same 
role of the entalpy in the phase transition process. 

We can assume that the function f{\/E 1 ) governs the first order phase transition from fermions (nucleons) to bosons 
(QD), capable of producing a compound system, which then proceeds to a saddle point, to scission and finally to 
fission. 

We can define the nuclear latent heat dL of the transition process as 



"-M^j ■ (23> 

From thermodynamic arguments and taking into account the transmutational potential defined in Eq.(9) and (11), 
the latent heat, at the thermodynamic equilibrium, can be written as 



T 



2 S N S D dn 2N _ 



N N N D dT 



(24) 



where S N and S D are, respectively, the entropies of the nucleons and the QD. 

From Eq.(24), we can realize that the important physical quantity is the variation of the transmutational potential 
with the nuclear temperature T (or with y/E~,); infact, this quantity modifies the balance of the entropy of the two 
phases of the system. n 2N _ tD (T) is an increasing function (except for 209 Bi, where rj 2N _ tD is nearly constant and the 
phase-transition is not relevant), therefore —dr\ 2N ^ D jdT decreases the latent heat, favoring the phase transition or, 
equivalently, increasing the photofission probability. 

IV. MICROSCOPIC INTERPRETATION OF THE TRANSMUTATIONAL POTENTIAL 

Let us now outline a microscopic interpretation of the transmutational potential. In the BCS theory, the ground 
state is composed by particles all paired with opposite spins: 
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\BCS>=Y[(u k +v k ala f _ k )\0> , (25) 

fc>0 

with the normalization condition: 

u 2 k + vl = l . (26) 

In Eq.(25), v\ and w 2 , represent, respectively, the probability that the state k is occupied by a pair of particles or not. 
We recall that the antisymmetrization is contained in the anticommutation properties of the fermionic operators and 
the product operator applies only over the k > states. The states k < refers to the conjugate states having the 
third component of the spin with opposite values. 

In our case, where QD pairs are contained into the nucleus, the state k can be empty or occupied by one non-paired 
nucleon or by a proton- neutron pair in a triple spin state. In analogy with Eq.(25), we can write the nuclear QD state 
at a given excitation energy in this form: 

\NQD>= Y[(u k + c k al + vl =1 ala f _ k )\0 > (27) 

fe>0 

where, now, —k indicates the state with third isospin component having the opposite value. 
The condition in Eq.(26) becomes: 

u\ + cl+v 2 k = l (28) 

where u\ is the probability that the state k is non-occupied; c 2 is the probability that the state k is occupied by one 
nucleon; v 2 is the probability that the state k is occupied by one QD pair. 

The three coefficients u k , c k , v k are functions of the excitation energy. With the above definitions we can write: 
and using the definition of the transmutational potential given in Eq. (9): 

2' 



V 2N ^ D (E 7 ) = k B T\og 



E 



(30) 



This potential contains all the physical information on the fermion-boson transmutation valid in the particular context 
of photofission of heavy nuclei, being fixed, at the moment, the region between 40 and 100 MeV. 
Note that if v k = (i.e. if inside the nucleus there arc no pairs contributing to photofission), the transmutational 
potential from two nucleons to a QD pair becomes equal to minus infinite, in agreement with the results one can 
deduce from Eq.(19) of Rcf. [10]: all the particles are fermions (nucleons). In conclusion, we remark that the state 
\NQD > of Eq.(27) does not represent, as the state \BCS > of Eq.(25), a redefinition of the vacuum state after 
having accomplished a Bogoliubov transformation; the state \NQD > represents an excitation state of the nucleus. 



V. CONCLUSION 



Among the several nucleon pairs that the correlations can arrange into nuclei, the quasi-deuterons are important in 
specific processes like, for instance, the photofission in the QD energy region. Infact, the empirical-phenomenological 
QD model describes accurately the behavior of the very different photofission cross sections of the nuclei 209 Bi, 232 Th 
and 238 U. 

This description can be theoretically justified on the basis of the statistical transmutation fermion (nucleon) - boson 
(QD) due to the transmutational potential studied in this work. The goal has been accomplished by considering, first, 
a semiclassical kinetic model which describes the distribution functions of nucleons (fermions) and quasi-deuterons 
(bosons) and the transmutation from one statistical family to the other one. These quantities are given in terms of 
the nucleon and QD chemical potentials, of the TP or the transmutation rate and are calculated in the framework 
of the QD model of photofission we introduced in the past and applied to the photofission of 209 Bi, 232 Th and 238 J7. 
We have found that the TP rj 2N ^ D (E 1 ) is related to a phase transition in the nuclear system. This has be shown by 
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means of a Clausius-Clapeyron equation satisfied by the photofission probability. 

Finally, we have indicated the existence of a relationship between statistical models, based on the semiclassical kinetic 
approach and many-body microscopic approaches like the superfluid nuclear model and have shown a microscopic 
interpretation of the transmutational potential. 
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Table Caption 

F 2 {E 1 ) for the three nuclei 209 Bi, 232 Th, 238 U. 
Figure Caption 

The transmutational potential as a function of the excitation energy for W9 Bi, 232 Th, 238 U (V„ = -48 MeV). 
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£ 7 [MeV] 


40 


50 


60 


70 


80 


90 


100 




0.0002 


0.0003 


0.0015 


0.004 


0.008 


0.012 


0.019 


•2-62 Th 


0.109 


0.142 


0.178 


0.221 


0.269 


0.324 


0.383 


238 jj 


0.183 


0.237 


0.296 


0.396 


0.428 


0.501 


0.578 
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